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ABSTRACT 

This paper is the first of a series dealing with the motion of flood waves and 
other waves of translation in open channels. The case treated is that of waves 
for which the forces of fluid friction are negligible with respect to the inertia 
and gravitational forces. The irrotational motion of a perfect liquid in a hori- 
zontal rectangular canal when the original surface is disturbed is investigated 
on the assumption that the horizontal velocity in a cross section is approximately 
uniform. The results are also applicable to motion in a canal of uniform slope 
containing water originally moving with a uniform velocity. Special emphasis 
is laid on disturbances which are propagated without change of form, and in 
these cases formulas are derived for the wave profile and velocity of propagation. 
Formulas are also derived which give the deformation, energy, motion of the 
center of gravity, and moment of instability of an arbitrary intumescence. Con- 
sideration is given to the maximum height of a wave of permanent form. For- 
mulas have been compared with the available experimental data. Of special 
interest is the comparison of the shape of the undulations composing the head of 
an initial sur6e with the characteristics of the cnoidal wave. 
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LIST OP SYMBOLS 

c velocity of propagation of long wave of negligible height and 

curvature. 
Ci velocity of wave approaching a discontinuity in a channel 

cross section. Also a constant of integration. 
c 2 velocity of transmitted wave after passing a discontinuity in 
a channel cross section. Also a constant of integration. 
Jacobian elliptic function, cosine amplitude, 
total energy of a wave per unit width of channel [MLT~ 2 ]. 
complete elliptic integral of the second kind. 
/ ( )> F ( )» w ^h or without subscripts — functional symbols, 
incomplete elliptic integral of the first kind, 
complete elliptic integral of the first kind, 
acceleration of gravity. 
g ( ), G ( ) functional symbol. 

height of main body of discharge wave (section V-17) ; vertical 
displacement of water surface. 
hi vertical displacement of surface as wave approaches discon- 
tinuity in channel cross section (section IV-5). 
h 2 vertical displacement of surface after wave passes discon- 
tinuity in channel cross section (section IV-5). 
hi maximum height of solitary wave, maximum depth of nega- 
tive wave, maximum height of initial discharge wave. 
hi, h 2 maximum height and depth, respectively, of cnoidal wave. 
hz=Hrj 3 . 

h a maximum height of initial discharge wave. 
h' average height of discharge wave. 
H undisturbed depth of liquid in channel. 
H ( ) functional symbol. 

k modulus of elliptic integral. 
I length of section of channel (section IV-4). 
l f m, n direction cosines. 

M 9 moment of instability [L], 
M' h Mi, M' 2) M x momenta per unit width [MT" 1 ]. 
p pressure [MLr 1 T~ 2 ]. 
p a atmospheric pressure. 
q speed of particle (section II— 5, 7, 8). 
q discharge per unit width in a cross section [L^T* 1 ] (section 

V-2). 
Q total volume per unit width of a solitary wave [L 2 ]. 
s length of arc. 
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sn (x, k) Jacobian elliptic function, sine amplitude. 
t time. 
u, v, w x component, y component, z component, respectively, of 
velocity. 
uq x component of velocity at bottom of channel, also first 
approximation to velocity in cross section. 
uo f *>o, wo velocity components prior to impulse (section II-6). 

tti, u 2 velocities of particles in a wave before and after passing a 
discontinuity in a channel cross section (section IV-5). 
tfi, u 2 , U\, U 2 local and mean velocities in two cross sections (section V-16). 
u\ velocity of particle at apex of wave (section V-19). 
U mean velocity of particles in cross section. 
Uq mean velocity in a section where h—Q. 

V volume [L 3 ] (section II-3); volume per unit width [L 2 ] (sec- 
tions IV-3, V-5). 
x Cartesian coordinate, specifically, parallel to channel in 

bottom plane. 
y Cartesian coordinate, specifically, lateral to channel in the 

horizontal plane. 
z Cartesian coordinate, specifically, drawn upward with origin 

at bottom of channel. 
a=3/H 3 (section V-10); a ratio (section V-14). 
7ii 72, 73 rates of angular dilatation [T" 1 ]. 
r circulation [LPT' 1 ]. 
ei, e 2 , €3 rates of linear dilatation [T -1 ]. 
e a small quantity. 
£, 77, £ components of vorticity [T" 1 ] (sections II— 1, 4, 5). 
f dimensionless x coordinate of free surface. 
77, £ z and x coordinates, respectively, of center of gravity of 
volume above undisturbed surface (sections V-5, 6, 12, 17). 
77 dimensionless z coordinate of free surface=A/H (sections V-8, 
9); Boussinesq number (section V-16). 
y', v", y'" — —m three roots of the differential equation of the permanent wave. 
77! maximum relative height of a solitary wave (section V-8); 
maximum relative height of a cnoidal wave (section V-9, 
18). 
% maximum relative depression of a cnoidal wave. 
9 ( ) functional symbol. 

X wave length of a cnoidal wave (sections V-9, 18). 
p density of liquid [ML -3 ]. 

<r partial volume per unit width of an intumescence [L 2 ]. 
S volume per unit width of an intumescence. 
<f> velocity potential [UT" 1 ]. 

4>o velocity potential prior to application of impulse (section 
II— 6) ; velocity potential, first approximation. 
^ ( ) functional symbol. 

co velocity of propagation of an element of volume of an in- 
tumescence, 
coo velocity of propagation of an element of wave height, negli- 
gible surface curvature, 
wi velocity of propagation of an element of wave height, appre- 
ciable surface curvature. 
ca e velocity of propagation of a cnoidal wave. 
a Q velocity of propagation of discharge. 
co impulsive pressure [ML _1 T~ 1 ]. 
ft gravitational potential [L 2 T -2 ]. 

I. INTRODUCTION 

The growing importance of predicting the occurrence of floods and 
the rate of travel and height of flood waves as they pass down rivers 
has led in recent years to a marked increase in the literature on this 
subject. However, in no language can there be found a comprehen- 
sive presentation of the mathematical treatment of the problem of 
such waves. Instead, various authors at different times and in dif- 
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ferent countries have treated separate aspects of the problem, so that 
the student of the subject must read many papers in several languages 
to obtain a satisfactory grasp of the problem. There is definite need 
for a systematic treatment of the results of these different methods 
of attack, presenting in coordinated form the results that have been 
attained up to the present time. 

It is for this reason that the National Bureau of Standards, at the 
suggestion of the United States Weather Bureau, has commenced the 
preparation of a series of papers dealing with the mathematical theory 
of flood waves and other waves of translation. It is not the purpose 
of these papers to furnish a practical method of predicting the rate 
of flood -wave travel and the rate of attenuation of the wave in an 
actual river channel with its complex flow conditions. Instead, the 
purpose is to furnish a sound mathematical theory on which attempts 
to solve the practical problems of flood prediction can be based. 

The present paper, the first of the series, deals with the motion of 
translation waves in channels of uniform, rectangular cross section 
when the forces of fluid friction are negligible with respect to the 
inertia and gravitational forces. Later papers in the series as now 
planned are in various stages of completion and will deal with the 
following topics: The effect of turbulence and channel slope and con- 
figuration on the motion of translation waves; the theory of quasi- 
permanent regime and the methods of prediction of flood waves ; and 
recent advances in the problem of the deformation of an intumescence. 

The different regimes of the flow of water in canals and natural 
watercourses may be characterized as uniform, varied, or unsteady. 

In the uniform regime the depth and mean velocity are invariable 
with respect to both time and space. This regime is attained in 
sufficiently long stretches of natural courses, provided the slope and 
discharge are constant, and the channel has a uniform cross section. 

In varied flow the depth and mean velocity are unchanging with 
respect to time, but not with respect to space. Examples are flow 
in the backwater region of a dam, flow in a canal connecting reser- 
voirs at different elevations, and flow in transition sections. 

On the other hand, the unsteady regime is characterized by a depth 
and velocity which vary with respect to both time and space. 
Gravity waves and capillary waves belong to this class of phenomena, 
the former being the more important of the two, considered from the 
point of view of applied science. Gravity waves fall into two classes — 
surface waves and translation waves. In surface waves, which occur in 
relatively deep waters, the motion decreases rapidly from the surface 
downward. The wave motion is sinusoidal, the individual particles 
periodically traverse closed paths, and the vertical accelerations are 
always appreciable. In translation waves, on the other hand, the 
motion of all the particles in a cross section is sensibly the same at 
any given instant. As a translation wave passes, the liquid particles 
are permanently displaced by an amount proportional to the volume 
of the wave. Waves of this type are frequently encountered in 
channels and rivers. Tidal waves are examples of the same phenom- 
enon in seas and oceans. 

Perhaps the simplest classification of translation waves is that 
obtained by considering the physical condition governing the motion. 
Two extreme cases are immediately apparent. First, we have the 
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case in which the effect of friction is negligible. In these waves the 
volume and the length are small; and the latter fact explains why the 
effect of friction can be neglected, so that the weight and inertia 
become the factors controlling the motion. If the curvature of the 
wave surface is negligible, the velocity of propagation of t he wave 
height, co , is given by the Lagrange velocity law, oo ^=^JgH y where g 
is the acceleration of gravity and H is the channel depth. On the 
other hand, if the maximum curvature and the height above the 
undisturbed level of the free surface are appreciable, then the velocity 
of propagation of the wave height «i equals co times a factor slightly 
different from unity, the exact value of the factor being a function 
of the height and the derivatives of the wave surface. We call such 
disturbances potential intumescences or irrotational translation waves. 
They can be generated by suddenly immersing a solid body in the 
liquid or suddenly adding a small mass of liquid; they may also be 
formed by the sudden motion of a body immersed in the liquid. In 
these cases the disturbed position of the surface is entirely above the 
original surface; that is, the waves have a positive elevation. The 
solitary wave is a familiar example. Waves can also be produced by 
abruptly removing from the channel a small quantity of water or an 
object already in the channel. These waves have a negative ele- 
vation. Surges produced in channels by an abrupt change of flow at 
either end also belong in this category. 

At the other extreme, we find the quasi-permanent regime of rivers, 
or what are commonly called flood waves. These are also translation 
waves, but they are waves in which friction is the all-important factor, 
the effects of weight and inertia being negligible. The waves are very 
long, which accounts for the predominant effect of friction. It is 
illuminating to compare flood waves and potential intumescences on 
the basis of their velocities of propagation. When discussing flood 
waves, it is convenient to utilize the idea of the velocity of propagation 
of discharge, <o tf . According to the Kleitz-Seddon law [1] l this velocity 
is given by dq/dH, the slope of the rating curve, q, being the discharge 
per unit width occurring at the normal stage, H. In a rectangular 
channel co ff =3i7/2, if the Chezy resistance law is used. A comparison 
of the two wave velocities gives co 2 lo) q 2 =(4:l9)((jH/U 2 ). Since in 
rivers the ratio gH/U 2 is a large number, the conclusion is that o? is 
many times larger than u q . This difference of velocities is significant. 

It should be understood that in reality there exists a continuous 
transition between the two extreme types of waves described above. 
There are other phenomena occupying an intermediate position, in 
which the effects of inertia, height, and friction are of the same order 
of magnitude. These phenomena are also of engineering interest, par- 
ticularly in connection with hydroelectric plants. When the turbines 
of such plants are suddenly shut down, a positive intumescence is 
generated in the supply canal and travels upstream, increasing in 
length as it goes. A negative intumescence is simultaneously pro- 
duced in the discharge canal, and this travels downstream, decreasing 
in length. These waves may be of considerable height (1 meter or so) 
and are attenuated by the forces of friction. Friction cannot be 
neglected, because of the long distances which the waves travel. 

» Figures in brackets indicate the literature references at the end of this paper. 
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Similarly, when the turbines are suddenly started, the same phenome- 
non occurs, but the signs of the intumescences are reversed. If the 
pressure and velocity variations caused by these waves are large, 
these effects must be considered in designing channel walls, bottom, 
and embankments. Waves like the above, except that they are 
generally propagated in still water, occur in canal reaches adjacent 
to locks when the latter are operated. The positive and negative 
waves that are generated are quite long, traverse large distances, and, 
although they undergo a deformation of form owing to the action of 
friction, the wave length remains the same [2], There are also a few 
special cases of flood waves occurring in mountainous regions, in 
which the inertia effects are important. Mass6 [3] has given the 
name of average intumescences to the phenomena described above. 

The present paper deals exclusively with potential intumescences. 
Since friction is negligible for these waves, and since they are generated 
in a liquid either initially at rest or in uniform motion, then by the 
Lagrange-Cauchy theorem the motion of the liquid particles is irro- 
tational. This fact justifies the mathematical approach selected for 
this paper. 

A brief presentation of the elements of irrotational theory is given 
in section II, primarily to emphasize the basic dynamical and kine- 
matical principles involved. Section III is a formulation of the basic 
differential equations whose solutions constitute the solution of the 
problem. 

The general theory is applied in sections IV and V, which deal with 
solutions of the equations to a first and to a second approximation. 
In the first approximation, the height and curvature of the waves are 
assumed negligible, and a complete solution is found in which the 
waves are propagated without change of form. In the treatment of 
the second approximation, which constitutes the main body of the 
paper, more general assumptions are made, leading to involved differ- 
ential equations which can be solved easily only for particular cases. 
The deformation of the wave profile during propagation and the motion 
of the center of gravity of the intumescence play an important role 
in this analysis. These points are discussed after the introduction 
of the concept of velocity of propagation of a volume element, origi- 
nally due to Boussinesq. There are three integrals invariant with 
respect to time associated with intumescences of appreciable curva- 
ture — the volume, the energy, and the moment of instability — which 
are discussed in considerable detail. 

The main emphasis is on the intumescence whose wave profile is 
propagated without change of form, and the particular solution is 
obtained for this case. If the liquid is undisturbed at infinity, we 
have the well-known solitary wave which has been investigated 
experimentally. If the liquid is in motion at infinity and this motion 
is indeterminate, we have the less familiar cnoidal waves. In the 
former case, there is an interesting principle discovered by Bous- 
sinesq — namely, for all intumescences of equal energy, the solitary 
wave is unique in possessing the minimum moment of instability. 
This theorem is significant in that it explains the frequent occurrence 
and ease of formation of the solitary wave. 

In discussing the velocity of propagation of surges, an attempt is 
made to determine the effect of the channel velocity distribution. 
The characteristics of the initial wave in a surge are not without 



pa«er?<m] Irrotational Translation Waves 53 

interest, and this wave is considered in the light of both theory and 
experiment. 

The concluding section touches on the problem of limiting wave 
heights and the breaking of waves. 

In the preparation of this paper, extensive use has been made of 
Lamb's comprehensive text [4] and of a classic paper by Boussinesq [5]. 
Whenever appropriate experimental data have been available, the 
authors have attempted to discuss the agreement between theory and 
experiment. The principal sources of data have been the well-known 
researches of Bazin [6] and Russell [7, 8] and a recent paper by 
Favre [2]. 

II. IRROTATIONAL MOTION 

1. DYNAMICAL EQUATIONS OF MOTION 

We consider an elementary parallelepiped 8x8y8z having its corner 
nearest the origin of the rectangular axes of reference at the point 
(x, y, z). Designating the component velocities at the point (#, y, z) 
by u, v, w, and the density of the liquid by p, we have for the kinetic 
reaction in the x direction of the liquid in the parallelepiped 

P^8x8y8z } (1) 

which is equal to the x component of all the forces acting on the liquid 
in the parallelepiped. If the pressure at the point (x, y } z) is p } we 
have for the x component of the force due to pressure, 

-^8x5y5z. (2) 

Denoting the potential of gravity at (x, y, z) by Q, we have for the x 
component of force due to the weight of liquid, 

- P ^8x8y8z. (3) 

Equating eq 1 to eq 2 plus eq 3, dividing by p8x8y8z } and applying the 
same reasoning in the remaining two directions, we obtain 

du__ 1 bp dfl 
dt~~ p bx bx' 



(4) 



dv__ 1 bp dO 

df~ p by by 

and dw_ 1 bp dO 

di~~ p bz bz 

To introduce the local variations of the velocity components, we 
consider u, v, w as functions of x, y, z, and t . With this understanding 
we write the identity 

dt~~ bx dt + by dt + bz dV bt' W 

where dx/dt } for example, represents the time variation of the position x 
of the liquid parallelepiped; that is, the velocity component u, and 
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bujbt represents the time variation of u at the fixed point {x, y, z). 
Interpreting the quantities dy/dt and dz/dt in the same manner, we 
write from eq 5: 

du bu , bu , bu , bu 



Similarly, 
and 



dt bt ' bx ' dy ' £>2 
ofe dtf , d# , dw . bv 

dw bw , bw , bw , dw 



(6) 



It should be noted that du/dt represents the acceleration of a particular 
element of the fluid followed in its motion relative to the axes x, y, z. 
The symbol d/dt has in general this significance. 

2. DISPLACEMENT AND DEFORMATION OF A LIQUID ELEMENT 

We consider the same elementary parallelepiped, 8x 8y dz f having 
its corner nearest the origin at the point {x f y, z). Designating the 
component velocities at some instant at the point (x, y } z) by u, v, w, 
we find for the relative velocities at the point (x-\-8x, y-\-8y, z+Bz) at 
the same instant, by Taylor's theorem, 



bu 
bx 



bu 
by 



sv^izfx+jSy+T 1 **' 



bu. 



and 



If we write 



bv* , bv^ . bv, 

Sv =b-x Sx+ by 8 y+b- z 8z > 

. bw. . bw. . bw. 
8w== te Sx+ b-y Sy+ te Sz ' 



(7) 



Bl = 


bu 

~ bx' 


bv 


bw 






t\" 


bw bv 
~b^ + bz' 72 ' 


bu bw 


bv 
~~bx 


bu 
+ by 


*= 


_bw 
"by" 


bv 
bz 


bu bw 
' bz bx' ^~ 


bv 
bx 


bu 
W 



and 



the relative velocities may be written 



and 



Bu=cifa+2(y£y+y£z) + 2^—^) > 



bv^^by+^a.bzAr^^+^bx—^z), 
dw=e^z+^{r2bx+yi8y)+-^(^y—Tidx)' 



(8) 



(9) 
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The following physical interpretation is now evident. The motion 
of an elementary parallepiped having the point {x, y, z) for its corner 
nearest the origin may be considered as made up of three parts. The 
first part represents a motion of translation of the whole element with 
the component velocities u, v, w. The second part represents the 
relative displacement of the particles within the parallelepiped and 
describes its deformation. The linear dilatations ei, e 2 , € 3 give the rate 
of increase of the edges of the element bx, 8y, 8z, respectively, during 
the motion. The angular dilatations y iy y 2y 73 give the rate of change 
of the angles of the edges of the parallepiped, taken two at a time. 
The third part represents a rotation of the whole element about an 
instantaneous axis having as components the angular velocities 
(1/2) £, (1/2) 77, (1/2) f. The vector (£, rj, f) is referred to as the vor- 
ticity of the liquid at the point (x, y, z). 

When, throughout a finite portion of a fluid mass, the component 
vorticities £, rj, f vanish, the motion of that portion is called irrotational. 

3. KINEMATICAL RELATIONS 

Since €1, e 2 , € 3 are the rates of dilation of the edges of a cubical ele- 
ment, it can be readily shown that the change, dV, of the volume of 
the liquid mass, AV= 8x by bz, due to the displacement of the particles 
per unit time is 

dV=(e l + e 2 + e 3 )V. (10) 

We also have the relation, 

Po V=(p Q +Ap)(V+dV), (11) 

derived from the principle of conservation of mass. p is the initial 
density of liquid and Ap the change in the density during the motion 
of the element. If the liquid is incompressible, Ap=0, and we obtain 
from eq 11, dV=0 f and from eq 10 

which is a kinematical relation that is satisfied in the region occupied 
by an incompressible liquid. Equation 12 is known as the equation oj 
continuity for an incompressible liquid. 

At the boundary of the liquid the equation of continuity is replaced 
by a different relation. We shall consider the relation in its most 
general form. Let F {x, y, z> t) = be the equation of the bounding 
surface. Let u, v y w be the velocity components of a liquid particle 
on the boundary, having the position (#, y, z) at time t. At time 
t-\-At the position of the particle is 

(x-\-uAt, y+vAt, z+wAt). 
Since the particle does not leave the boundary [4, p. 7], we must have, 
because of the equation of the boundary, 

F(x+uAt, y+vAt } z+wAt, t+At)=0. 
When the latter expression is expanded by means of Taylor's theorem, 
we obtain the relation 

dF bF . Z>F t bF. bF . , 1Q . 

W=M +u lte + v bTj +w b-z =Q > (13) 

which is the general boundary condition. 
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Two special forms are to be noted. At a fixed boundary 

lu+mv+nw—O f (14) 

where I, m, n represent the direction cosines of the normal to the 
boundary. At a surface of discontinuity 



l{u x — u 2 )+m(vi— v 2 )-\-n(wi— w 2 )=0, 



(15) 



where the suffixes are used to distinguish the velocities at the two 
sides of the surface. 

4. IRROTATIONAL MOTION 

In a large and important class of hydrodynamical phenomena, the 
vorticities £, tj, f vanish — that is: 



k- 



bw 
'by' 



bv 
bz s 



=0, 



and 



bu bw 



r= 



bv 
bx 



dy 



(16) 



As a result of this restriction, the expressions for the velocity com- 
ponents, u, v, w, are very much simplified. :* ->, 

We consider now the function 0, defined by f the' differential 
equation 

—d<t>=udx+vdy+wdz. (17) 

According to a theorem in differential equations, the necessary and 
sufficient conditions for d<f> to be an exact differential, that is, for 
d<j> to be expressed as 

are the relations in eq 16. Comparing eq 18 with eq 17, we write for 
the velocity components of 'a liquid particle of zero vorticity: 

50 

bx* 



u=— - 



50 



and 



w=* 



50 
bz 



(19) 



The function is called a velocity potential from its analogy with 
potentials arising in some other branches of physics. As the velocity 
potential exists if and only if the flow is irrotational, the latter flow is 
also called potential flow. 

?•■ If the liquid, in addition to being in irrotational motion, is also in- 
compressible, [then the velocity potential satisfies Laplace's 
equation, 



5 2 5 2 5 2 0_ 
bx 2+ by 2+ bz 2 



0. 



(20) 
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This can be seen by substituting the expressions for the velocities 
given in eq 19 in the equation of continuity, eq 12. Thus, for 
the conditions of flow assumed here, (the irrotational motion of an 
incompressible liquid) the complete determination of the velocities 
at all points reduces to the selection of the harmonic function <j> 
satisfying the boundary conditions. 

5. INVARIANCE OF THE CIRCULATION IN A MOVING CIRCUIT 

The following theorem on irrotational motion is of great importance; 
If the motion of any portion of a perfect fluid is irrotational at any one 
moment, it will continue to be irrotational, provided the external forces 
have a potential and the density of the fluid be either constant or a function 
of pressure only. 

One of the methods for establishing the theorem is to consider the 
circulation around a closed carve of particles and inquire into the 
variation of the circulation as the particles of the curve move into new 
positions. By the circulation around the circuit, G, is meant the line 
integral at any time, t, 

(udx+vdy+wdx), (21) 



i 



of the velocity vector^ taken around C, a closed path of particles in the 
liquid. The circulation of the circuit of the same particles at time 
t+dt is 

where 

^= <j) ^(udx+vdy+wdz). (22) 

To evaluate dYjdt we take an element of circuit Oof length 8s and con- 
sider the time-rate of variation of the quantity uhx as the particles on 
the element move into new positions. Differentiating the product 
uhx, we obtain 

£(tto) =£**+«£ to. (23) 

The first term in the right-hand member contains the x component 
of the acceleration of a particle on ds, given by the first of the following 
dynamical equations of motion 



and 



du_ 1 dp 
dt~ pi>x 


dx f 


dv__ 1 c)p_ 
dt~~ p by 


~dy 


dw__ 1 dp__ 
dt~ p dz 


da 

bz t 



(24) 
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Here p is the pressure, p the density of liquid, and ti the potential 
of the external forces at the point under consideration. The second 
term is readily transformed into 

u j. (8x)=u8 -j=uZu, (25) 

since dxjdt is the velocity u. Thus, 



and similarly 



and 



d . . . 1 dp dfi 



(26) 



Introducing eq 26 into eq 22, we obtain 

where ([ 2 =u 2 +v 2J rw 2 . Since 12 and q 2 are single-valued, continuous 
functions of x, y, z, the last two integrals in the right-hand member 
of the equation vanish. If p is a single- valued, nonvanishing, con- 
tinuous function of p only, or a constant, then the first integral 
vanishes also. Hence, under these circumstances 

f=0; (28) 

that is, the circulation around a circuit moving with the liquid re- 
mains unaltered. 

By using Stokes' theorem, the circulation, T, around circuit C, 
may be related to the vorticity of the liquid particles on a surface S 
bounded by 0. Stokes' theorem states that the line integral of any 
vector around a closed curve is equivalent to the flux of the curl of 
the vector through any surface bounded by the closed curve. Using 
this theorem, we now write 

T = f s (li:+mr J +nt)dS, (29) 

where Z, m, n are the direction cosines of the normal drawn inward 
from the surface. 

It now follows from eq 28 and 29 that the property of irrotation- 
ality of a portion of the fluid is carried with the fluid. For, if initially 
£, 77, f vanish in a region containing a fluid, then the circulation for any 
circuit in the same region vanishes. This is true at all subsequent 
times in the region containing the same particles, by eq 28 above, and 
the vorticities consequently remain equal to zero. 
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The restrictions under which the main theorem, expressed by eq 28, 
has been proven must be remembered. It is assumed that the 
externa] body forces have a potential, and that the density is either 
constant or a function of pressure only. The latter condition is 
violated, for example, in the case of convection currents generated by 
a temperature gradient. Since in the dynamical eq 24, which were 
employed in the proof, the viscosity terms are missing, another 
restriction is the condition that the effects of viscous forces are 
negligible in comparison with inertial forces. 

6. IMPULSIVE GENERATION OF MOTION 

When the boundary conditions of a mass of liquid are suddenly 
changed, a sudden alteration in the motion may take place. The 
forces introduced are of the nature of impulsive pressures, and, during 
the time interval that the impulses last, the effects of all the other 
forces are negligible. 

Let p be the density, u , v , w the components of velocity imme- 
diately before, u', v' , w r those immediately after the application of the 
impulse, Z> the impulsive pressure at the point (x, y, z). The change 
of momentum parallel to x of an elementary rectangular parallelepiped 
having edges 8x, 5y, 8z, and with its corner nearest the origin at 
(x, y, z), is p5x5ydz(u'--~'Uo); the resultant impulsive pressure in the 

same direction is —-^—8x8y8z- Hence, 

pbxhybz (V — u ) = — j—bxbybz, 
or 



Similarly 
and 



u f - 


-u =- 


_i as 

p da: 


v'- 


- v Q =- 


_1 dw 

p by' 


w'- 


-vj =- 


_1 dw 

p dz , 



(30) 



In these relations we have a physical interpretation of the velocity 
potential, <£. For, consider the case where Uo, v , w are zero. A 
state of motion for which a velocity potential exists could be produced 
from rest by the application of a properly chosen system of impulsive 
pressures. If <j> be the velocity potential, the desired system of impul- 
sive pressures is o?=p0+(7. The arbitrary constant in this expression 
shows that a uniform pressure throughout a liquid mass has no effect 
on the motion. If the initial state of motion be uniform, that is, 
if for the initial state of motion the liquid possesses the velocity 
potential 9 <t>o=UoJc-\-VQy+w z , then the new state of motion possessing 
the velocity potential <f> may be produced from the initial motion by 
the application of the impulsive pressures w=p(#— 0o) + C 



~1 
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7. PRESSURE EQUATION 

In a case of irrotational motion the velocity potential, <j>, determines 
the velocity components u } v, w, but the hydrodynamical solution of 
the flow is not complete as long as the pressures are unknown. The 
determination of the pressures is made by means of the pressure 
equation, which is a dynamic relation to be derived presently. 

Since the motion is irrotational, the velocities u> v, w in the dynamical 
equations, 

du__ 1 bp 512 
dt ~~ p bx bx 



and 



dv 
dt" 


p <>y~ 


50 


dw 
dt ~~ 


p bz 


512 

""5?. 



(24) 



where the operator d/dt is 

d b . d . 5 . d 

may be expressed in terms of the velocity potential, <£. 
equations are transformed into 



5 2 <£ , 1 d 



512 



bxbt 

5 2 i<*> 



^uw^l 5 p 
2 bx bx p bx 



h il«. 



and 



dydr 2 5i/ 
5 2 <*> . 1 5 



r- 



50 



1 5g 

p &y : 



512 1 5^ 



5^5* ' 2 5s 



9 7\JL >v« 



5s 



5zJ 



(31) 
The dynamical 



(32) 



Multiplying the first by dx, the second by dy, the third by dz, adding, 
and then integrating the sum, we obtain the integral 



s 



dp_ d<ft 



L_Q_(i/2)^+y(<), 



(33) 



where f(t) is an arbitrary function of time. This equation makes it 
possible to determine the pressure, once the velocity potential is 
known; it is called the pressure equation. 
For an incompressible liquid, 



MH-k+^ 



(34) 



To determine the arbitrary function /(<) we must know p and <t> as 
functions of t at some definite point in the region occupied by the 
liquid. Generally this point is one where i>(j>jbt vanishes, q is a con- 
stant, and p is the atmospheric pressure p„. 
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8. SOLUTION OF PROBLEMS IN IRROTATIONAL MOTION 

The mathematical procedure for solving a problem of the irrota- 
tional motion of an incompressible liquid bounded in part by fixed 
boundaries and in part by a free surface consists in finding a har- 
monic function, <j>(x,y,z,t)y which satisfies Laplace's equation (eq 20) 
in the region occupied by the liquid, and the following surface condi- 
tions. 

At the free surface, if it exists, there are two conditions to be satis- 
fied. If F(x,y,z,t) is the equation of the surface, we have 

bF, bF. bF, ZF n , 1Q , 

where u, v, w are the velocities at the surface. The second condition 
states that the pressure is atmospheric at the free surface; we then 
have from eq 34 

2Wp=^-n-^z 2 +/(<). (35) 



At the fixed boundaries, 



|-0, (36) 



n being the normal drawn to the boundary. 

The function, <£, which satisfies eq 13, 20, 35, and 36 is unique, 
except for an 'additive constant. For references to proofs see [4, 
p. 59]. 

III. FORMULATION OF THE WAVE PROBLEM— WAVES IN 
WATER OF SMALL DEPTH 

We consider a rectangular horizontal channel of uniform cross- 
section and infinite length, containing a perfect liquid initially at rest. 
We assume that some external agency acting for a short time produces 
a deformation of the free surface of the liquid in the channel. If the 
action of the external agency is such that the associated forces exerted 
on a portion of the water in the channel are essentially in the direction 
parallel to the direction of the channel and are constant throughout 
the whole depth of the liquid, the particles in a normal section under 
the deformed surface will be moving with a constant velocity. Under 
these conditions the deformed free surface constitutes a translation 
wave. 

The liquid of the channel being initially at rest, we infer from the 
theorems previously explained that the resulting motion of the liquid 
under the deformed free surface is irrotational and possesses a velocity 
potential, <£. Let the axis of x be drawn in the bottom plane of the 
channel parallel to the channel axis; let the axis of z be drawn upwards. 
Denoting the velocity components in the direction of x and z by u 
and w, respectively, we have 



62 Journal of Research of the National Bureau of Standards [Vol u 



and 



— §*, (37) 



—f z > < 38 > 



where is a function of x, z, and t and satisfies Laplace's equation 

The pressure equation (35) reduces to 

2=^-ffg-|(«»+W»)+F(0. (40) 

where # is the intensity of gravity. 

Let the z coordinate of the free surface be denoted by h+H, where 
H is the initial undisturbed depth of water in the channel, and h is the 
elevation (or depression) of the wave. 

To evaluate F(t), let us consider the point at infinity on the free 
surface; at such a point d<l>/dt=u=w=h=0, z=H, and p=p a . Sub- 
stituting in eq 40 there results 

p 

At the surface, since p—p a , we obtain the dynamical surface condition 
from eq 40, neglecting the effect of surface tension: 

||— flrA— |(t*»+to») = 0, z=H+h- (41) 

Although the last equation has been derived for the case where the 
liquid at infinity is at rest, it is equally applicable to the case where 
u 2 -\-w 2 and b^fbt vanish at the points on the surface where h=0, as, 
for example, in an intumescence consisting of an infinite number of 
undulations of the same dimensions moving without change of form. 

We obtain the kinematic surface condition, that a fluid particle of 
the free surface always remains on the free surface, by substituting 
F(x,y,z,t) — z-—H—h=Q in eq 13, obtaining 

At the bottom surface of the channel the velocity component normal 
to the bottom vanishes, that is, at 2=0, 

g-0. (43, 

In view of this last condition, we adopt for the potential <£ the power 
series 

*=SAiS", *i=0, (44) 

n=0 
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where <j> n is a function of x and t only. Substituting this expression 
for $ in eq 39, and equating the coefficients of z n to zero, term by term, 
we arrive at the convergent series 

. z 2 d 2 tf> , £ 4 d 4 4>o 2 6 d 6 4> , fA§Ci 

Differentiating with respect to 2, 

d<ft b 2 <t> 3 3 d 4 <fo z 5 d 6 <ft 

dz~ J dz 2 +3! dx 4 5! dz 6 + " ' * * l4b; 

The function O , as stated above, depends only on x and t. If we 
denote the velocity of particles of water at the bottom of the channel 
at x by u , then 

*>=(*>).- f - ^<fe=(^)«+J^t^, (47) 

where (</> )» is the value of <£ at x=co . Since the liquid in the channel 
at infinity is at rest, we may take the value of (</>o)« as zero: then, 

$0= J u dx. (48) 

In the above formulation, h and <£ are functions specifying a trans- 
lation wave. They can be determined by integrating the two simul- 
taneous differential eq 41 and 42 representing the dynamical and 
kinematical surface conditions. 

IV. LONG WAVES OF NEGLIGIBLE HEIGHT AND 
CURVATURE 

1. VELOCITY OF PROPAGATION 

If u 2 and w 2 are negligible in comparison with gh (assumptions I, 
II) and udh/px in comparison with dh/dt (assumption III), then the 
surface conditions given by eq 41 and 42 reduce to 



and 



0»-D=O, (*» 



If+lH- <w 



Also, if in the expansions of <f> and d^/ds given in eq 45 and 46, only 
the first term is retained (assumption IV), we have 

*=$o, (51) 

and 

IS— I* < 52 > 

198881—40 5 
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Here, then, we have the basis for a first approximation to a solution of 
a wave configuration. The horizontal component of the velocity in 
a normal section is constant and equals the bottom velocity u ; 
the vertical component varies linearly with z. 

Substituting eq 51 and 52 in eq 49 and 50, we obtain 



and 



Sl-^-0, (53) 



assuming that h is small in comparison with II and hence can be 
neglected (assumption V). Eliminating <£ from eq 53 and 54, we 
obtain the differential equation 

5%_ 2 d 2 h ( . 

dt 2 ~ C dx 2 ' [bb) 
where 

c=<y/gIL (56) 

Eliminating h similarly and writing 

«.--!*. (57) 

there results for Uq the relation 

SMS?- < 58 > 

Equations 55 and 58, which are well known, possess complete 
solutions, which also satisfy eq 53 and 54: 

^F 1 (x-ct)+F 2 (x+ct) 9 (59) 

and 

/ ^=F 1 (x-ct)-F 2 (x+ct) } (60) 

c 

where Fi and F 2 are two arbitrary functions. 

These results have a very simple geometrical interpretation. 
Consider the motion represented by the function Fi in eq 60. When 
x and t are increased by ct and r, respectively, the value of the func- 
tion is not changed; the disturbance h which existed at the place x 
and at time t has moved during the time interval r to the place 
aj+cr. Therefore, the function F\ represents a wave traveling in the 
positive z-direction with out change of form at a constant velocity 
of propagation 2 c=^/gH. 

Similarly, F 2 represents another wave propagated in the opposite 
direction with the same velocity. Any motion whatever subject to 
the conditions assumed above for the approximate solution will be 
made up of waves of these types. F t and F 2 are determinate when 

* The question of the meaning of velocity of propagation when the wave does not retain its form will be 
discussed later in section V-l. 
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h and u are known for every value of t, and at a single point x=x x \ 
or when h and u are known for every value of x and for a single 
time t=ti. 

If h and u Q represent the surface elevation and velocity of particles 
in either the progressive or regressive wave (Ft or 2^=0) it is seen 
from eq 59 and 60 that 

Uo=±Cjj7 (61) 

where the plus sign appKes to a wave moving in the positive z-direction 
(i^ 2 =0) and the minus sign to a wave moving in the negative ^-direc- 
tion CFi-0). 

We shall now examine the restrictions assumed in obtaining our 
approximate solution: 

u 2 <zgh, (I) 

v?<gh, (II) 

and 

dh ^i>h /ttt\ 

The two series converge rapidly ; that is, 

z 2 £> 2 *o 
and 

ff.<l. (V) 

It can be demonstrated that the following physical restrictions lead 
to the above mathematical assumptions: 

(A) The ratio of the wave amplitude to the depth of the channel 
is small: (h/H)<$ 1. 

(B) The slope of the wave is small: (i>h/dx)<^l. 

(C) The radius of curvature of the waves is large with respect 
to H 2 l2h. 

2. ENERGY OF AN INTUMESCENCE 

A single wave is called an intumescence. The potential energy 
per unit width of an intumescence, at time t, due to the elevation or 
depression of the fluid surface above or below the level of the previ- 
ously undisturbed liquid in a channel is 



2! *> <*» < IV ) 



&■ f * h*dx t (62) 



where x t and ck 2 denote the extremities of the intumescence. The 
kinetic energy per unit width of the liquid under the free surface of 
the intumescence is 



*~7 1 v%dx 



(63) 
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Substituting for Uq its value in eq 61, eq 63 becomes 



29 
2 



f X% h 2 dx. (64) 

JXi 



Comparing eq 64 and 62, it is seen that the potential and kinetic 
energy of an intumescence are equal and that the total energy is 

gp ( X% h 2 dx. (65) 

The integral in eq 65 without the multiplier gp will be of special 
interest, and will be subsequently referred to as the energy, E } of 
an intumescence. 

3. DISPLACEMENT OF THE PARTICLES 

The liquid particles in the cross section at x, before the arrival of 
an intumescence, are at rest. At the moment when the intumescence 
is passing, the particles in the section are moving with the velocity u , 
to a first approximation. Thus, the volume of flow through the 
section during the passage of the intumescence is 

V= p u Hdt. (66) 

Since 

c dt=dx, (67) 

we have 

V =£ ^Ldx. (68 

Substituting for u from eq 61, we obtain 

V= ( Xi hdx. (69) 

JXi 

The volume of flow through a cross section equals the volume of the 
intumescence. Hence, during the passage of an intumescence, the 
particles in a cross section are moved forward a distance 8x, the mag- 
nitude of which is given by the ratio of the volume of the intumescence 
to the area of the channel cross section. 

4. EFFECT OF AN ARBITRARY INITIAL DISTURBANCE 

If for £=0 the velocity of the liquid particles in the channel and the 
form of the disturbed surface are known: 



c 
and 



*-»(*>, -•<*< + •, (70) 



ji~m, -co<s<+«, (7i) 
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then the functions F\ and F 2 occurring in eq 59 and 60 are 

(72) 



Fi(as)=|[*(*)+*(*)], 



and 

TO=f [-*(*)+*(*)] 

To determine the form of the wave profile, h/n—\l/(x,ti) for any 
subsequent time, t=t 1} the curve representing Fi(x) is displaced in the 
positive ^-direction by a distance ct u the curve representing F 2 (x) by 
the same distance in the opposite direction, and the ordinates of the 
two displaced curves are added algebraically. 

Consider disturbances confined initially to a finite segment of the 
channel of length I. Resolution of the initial disturbance into F\ and 
F 2 indicates that after the lapse of the time interval r=c/2, two waves 
of equal length I will be traveling in opposite directions. In general, 
resolved waves will have unequal heights at corresponding points. 
If, however, in the initial state, u—0, that is, 0(z) = O, then JF\(oO = 
F 2 (x), and the two component waves have the same shape. If, on 
the other hand, u= ±ch/H initially, then either F x or F 2 is zero, and 
the result is a single wave traveling in one direction. 

5. REFLECTION AND TRANSMISSION OF WAVES 

A wave is reflected when it reaches a barrier in a channel. We take 
the origin, x=0, at the barrier and note that u=0 at the barrier. 
Utilizing this fact in the general eq 60, we obtain the condition that 

F l (-ct)=F 2 {+ct) (73) 

for the incident and reflected waves. 

Applying this last condition to the general eq 59 and 60, we obtain 

~=F 2 (ct-x)+F 2 (ct+x), (74) 

and 

■^=F 2 (ct-x)-F 2 (ct+x) > (75) 

which are the equations for reflected waves or for identical waves 
traveling in opposite directions. Here, if F 2 (ct-\-x) is taken as the 
incident wave, F%(ct—z) will be the reflected wave. At x=x h and 
t=t u the elevation or depression of the wave is 

-^=F 2 (ct 1 +x 1 ). (76) 

The corresponding elevation or depression of the reflected wave will 
reach the point x=Xi at the time t=ti + 2xi/c 9 and the value of the 
elevation or depression is 

^=F 2 (c{t l + 2x 1 lc}-x 1 )=F 2 (ct 1 +x 1 ). (77) 
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Hence, elevations and depressions are reflected without change. It 
can be shown in a similar manner that corresponding velocities are 
reversed with their magnitudes unchanged. 

When a wave reaches a point in a channel where there is an abrupt 
change in the cross section of the channel, either in width or in depth, 
a part of the approaching wave is transmitted and a part is reflected. 
Taking the origin at the point of the abrupt change of section, we 
may write from the general eq 59 and 60 for the side from which the 
wave approaches, moving in the positive ^-direction, 

j^Fix-cd+ficJ+x), (78) 

and 

^Fix-crf-jfat+z), (79) 

in which the first term in the right member of each equation applies 
to the approaching wave, and the second term to the reflected wave. 
For the transmitted wave we have 

jf-gix-cj), (80) 

and 

~=g(x-c 2 t). (81) 

c 2 

In the above equations fli and H 2 are the depths of the channel 
on the two sides of x=0, and c x and c 2 the velocities of propagation 
for the corresponding portions of the channel. At x=0, hi=h 2 , since 
the surface is continuous, and therefore, 

Hm-Cit)+f(c x t)]=H 2 g{--c 2 t). (82) 

Also, the conservation of matter requires that BiHiUi—B 2 H 2 u 2 , where 
B x and B 2 are the widths of the channel on the two sides of x=0. This 
condition gives 

B 1 H 1 c 1 [F(-c l t)-f(c 1 t)]=B 2 H 2 c 2 g(-c 2 t). (83) 

We deduce from eq 83 and 84 the ratio of the elevations in the cor- 
responding parts of the reflected and incident waves and of the trans- 
mitted and incident waves. These are 



f ^ Bid—B^ 
F B x c x +B 2 c 2 

and 

£_ 2B l c 1 H 1 
F~(B l e l +B a fi 2 )H 2 9 



(84) 



(85) 



respectively. 

A wave traveling in a channel connected to a body of water in a 
basin of infinite depth and width may be considered as a special case. 
For this case it is evident by inspection that g/F=0 &ndf/F= — l. 
Hence the incident wave is completely reflected at the juncture, and 
the elevations are reflected as depressions. 
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V. WAVES OF APPRECIABLE HEIGHT AND CURVATURE 
1. TWO DEFINITIONS OF VELOCITY OF PROPAGATION 

For long waves of infinitesimal height, the concept of velocity of 
propagation entails no ambiguity, because such waves are propagated 
without change of form. However, when a traveling wave is under- 
going deformation, the meaning of velocity of propagation must be 
clearly defined. We may $dopt either of two distinct definitions. 
The first definition will be developed on a purely formal basis. In a 
given wave system we consider the quantities dh/dt and bh/bx for an 
arbitrary x and t, and write the following relation from eq 13: 

§+«i§=0,A=ACM). (86) 

In effect, this relation defines the quantity «i, which we may call the 
velocity of propagation of a wave-height element. When the traveling 
wave is undergoing a change of form at a given instant, t y co t is a 
function of h. 

A second definition, which is due to Boussinesq, considers the vol- 
ume of an intumescence. The volume of an intumescence is the 
quantity of liquid in the entire intumescence between the free sur- 
face and the undisturbed primitive level of the liquid in channel. We 
have, therefore 

" hdx f (87) 



Jxo 



where x is the rear of the intumescence, assumed moving in the 
positive ^-direction. The intumescence may be limited posteriorly; 
that is, for all values of x^x Q , h is negligible, and in this case 2 has a 
finite value. On the other hand, the intumescence may be unlimited; 
in this case we consider the volume, 2, included between the head 
and an abscissa, x , considerably removed from the head, and moving 
in such a way that 2 remains constant. 

We now consider a partial volume, <r, of 2, between the head of the 
intumescence and an abscissa, x: 



-r 



hdx,z£xo. (88) 



We assume that a plane originally at x moves in such a way that it 
always has in front of it the same partial volume a of the intumescence. 
We refer to the velocity at which this plane moves as the velocity of 
propagation of an element of volume of an intumescence, and denote 
this velocity by «. 

The relation between w and h may be obtained from eq 88 by 
differentiating the right-hand member of the equation with respect 
to t and equating to zero. However, it may be more instructive to 
obtain the desired relation in the following manner (see fig. 1). The 
curves ahc and a'b'c' represent the free surface of the moving intu- 
mescence at the times i and t+At, respectively. A plane at x having 
in front of it, at time t, the partial volume, <r, of the intumescence, 
moves to the position z+uAt at time t+At and still has the same par- 
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X+CJAt *"** 

Figure 1. — Diagram illustrating the velocity of propagation of a volume element. 



tial volume, <r, in front of it. The areas A and B are equal to each 
other. Consequently, 

ha>At= j ^Atdx. (89) 

Dividing by At, 



h:***** 



(90) 



and then differentiating with respect to x, we obtain 

which is the relation connecting o> with h. 

Eliminating dA/d£ between eq 86 and 90, we find that the difference 
between the two velocities of propagation is 



jdoo Ibh 



(91) 

If a wave is moving without change of form, o> is constant, and the two 
definitions give identical values for the two velocities of propagation, 
a fact that is self-evident. 

2. MEAN VELOCITY OF PARTICLES 

If U and g, respectively, denote the mean velocity of the particles 
in, and the discharge through, a cross section at x, then 



q=U(H+h) 



(92) 



Consider the inflow and outflow in a time interval, At } for the vol- 
ume between two cross sections separated by a distance, Ax; Ax 
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and At having the same order of magnitude. The mean discharge 
at x in the time At is 

. At dq: 
* + 2 *' 

omitting infinitesimals in At of the second and higher orders. The 
volume of the inflow in the time At is then 

., , (At) 2 bq 

<L At +—-it' 

Similarly, the mean discharge at x + Ax in the time At is 

2+ dx + 2 dt 

omitting infinitesimals in Ax and At of the second and higher orders. 
The volume of the outflow at x-\-Ax in the time At is 

qAt + £ A X At + - Y .^, 

and the accumulation of liquid in the volume is 

-^AxAt. 
OX 

The increase of volume of the liquid is also 

oh A A, 

bT AxAt > 

omitting infinitesimals of the same orders as before. 
Therefore, 

dg_ dh 

dz~~ d* 9 (93) 



or from eq 92 



|(l7Cff+tt]— jf- 



Substituting from eq 90 

|[tf(ff+A)]-|(M> (94) 

wliich when integrated becomes 

U(H+h)=ha+C. (95) 

Since the liquid is assumed to be undisturbed at infinity, the constant 
of integration vanishes, and the equation for the mean velocity is 

TT ho* 

U -H+h' (96) 
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3. EFFECT OF CURVATURE AND HEIGHT ON VELOCITY 
OF PROPAGATION OF A VOLUME ELEMENT 

When the ratio h/H and the curvature of the wave profile are appre- 
ciable, the first approximation is no longer satisfactory, and in such 
cases recourse to a second approximation becomes necessary. The 
method which will now^ be presented is due essentially to Boussinesq. 

In the surface conditions, eq 41 and 42, we neglect w 2 and replace u 
by the value of u in the first approximation, eq 61, replacing c by its 
value in eq 56: 

w=J|a. (97) 

These conditions now become 

flfr-H+g^-O, 08) 



and 

d<£ , bh , I g ^dh 



Retaining the first two terms in the expressions for <f> and d#/d2, eq 45 
and 46, and neglecting h in the terms having the smaller values, we 
obtain 



0=^0 — 2" "5^2-' ( 10 °) 



and 



We differentiate eq 100 with respect to t, replacing d<j&o/E>£ in the second 
term of the right-hand member by its value in the first approximation 
given by eq 53, that is, by gh. In this manner we obtain 

We next consider eq 101. In the right-hand member there are two 
small terms. In these terms we replace d<t> /dx by its first approxima- 
tion from eq 57 and 61, that is, by ~-W jM'- ^ n ^ n * s mai * n er, re- 
stricting ourselves to waves propagated in the positive ^-direction, 
we obtain 

Substituting from eq 102 and 103, we obtain from eq 98 and 99: 
and 

K^lg-f^o. (.05, 
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To eliminate <t> between these, we differentiate eq 104 twice with 
respect to x and multiply through by H, and differentiate eq 105 with 
respect to t. In place of actually differentiating with respect t o t, 
however^ we differentiate with respect to x and multiply 3 by —->]gH. 
Subtracting eq 105 from eq 104, there results: 

which is the fundamental wave equation to the second approxima- 
tion, in which the effects of the wave height and of curvature of wave 
profile are included. 

Although this equation has been derived for the case where the 
liquid in the channel is undisturbed at infinity, it remains v alid for the 
case where the liquid at infinity is in motion, provided u=^g/H h, to a 
first approximation (eq 97), 

Equation 106 becomes integrable once with respect to x when bh/bt 

is introduced from eq 90. Performing the integration with respect to x 

and remembering that h and its derivatives vanish for x= » t we 

obtain 

5,,, rjbh x rjb/M* , H 2 b 2 h\ n , 1AI7N 

Wt ( h „) + gH^+gH^ m+T ^ 2 J=0, (107) 

which permits the determination of oo, the velocity of propagation of a 
wave- volume element, in terms of h and d 2 h/bx 2 . 
We adopt as an expression for ca the product 

co=V^(l+^) J (108) 

where € is a small quantity. Multiplying the two members by h and 
differentiating with respect to t, we obtain, remembering that e is a 
small quantity and applying the same reasoning used in differ- 
entiating eq 105, 

Introducing bhfdt from eq 90, thisjuay be written, 

Jj<M=-V^(M-0ff|^*); d09) 

or, in view of eq~108, 

J(M = -<7tf^[(l+2e)A]. (110) 



8 The justification for replacing the operator b/dt by — *JgH bfdx can be seen from the following: Assuming 
that the wave retains its form, that is, that dFldt=0, where F is a function of (x—ct), we have the relation 



from which it is obvious that 



dF b_F bF 
dt mm bt+ c bx' m[i ' 



bt~ c bx M bx 



While the present discussion relates to waves that may change their form, the above relation can be applied 
to small terms in the wave equation, to a first approximation. 
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The substitution of eq 110 in eq 107 makes the latter expression 
immediately integrable with respect to x. Integrating, putting the 
constant of integration equal to zero, since h and its derivatives 
vanish at x= co, and then dividing by h y we find 

2e ~2H + Ud? (111) 

Hence 

which is the expression for the velocity of propagation of a volume 
element of an intumescence moving in still water. 

We next derive the expressions for u 0) the horizontal velocity com- 
ponent at the bottom, and u, w y the velocity components at a point 
having the ordinate z. To determine u we consider eq 105, replacing 
dh/dt by ~5(Aco)/dx from eq 90, and then integrating with respect to 
x, the constant of integration being equated to zero. In this manner 
we obtain, since -w =-—d<£o/2>£, 

•y^r s+Tw (113) 

From eq 45, after differentiating with respect to x, and from eq 46, 
retaining only the first two terms of the right members, we obtain: 

z 2 i> 2 Uo /11yl . 



u -- 



and 



bu , 2 3 d 3 u /t . K v 

tt =-*fe+6"5?-' (115) 



which give, after substituting for u from eq 113, and neglecting quan- 
tities of small order: 



u 



and 



w= 



^^H-m>+\^-m)-w\ J (116) 

,— T/l h \dh.l/rj z 2 \b*hl , 11in 

-^V^(^-2^^5+^^^l (117) 

4. DEFORMATION OF THE WAVE PROFILE 

The expression for o>, eq 112, and the relation in eq 90 permit the 
evaluation of the change in h occurring during a short interval of time 
to a sufficient approximation. Denoting -jgH by a? , eq 90 may be 
written as 

The integral of this expression is 

A=/(a:-«b0-*|^(«-M»)]f (11») 
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when t is small. This may be proved by differentiating eq 119 and 
substituting in eq 118. In the process of differentiating the second 
term in the right member, the relation b/dt=co i>/dx is used, this being 
valid if t is small. Eq 119 enables us to trace the wave profile curve 
as it would be for time ti+t, if it is given for time ft. 

We next take an element of the intumescence, do, having the 
height h and base dx and follow its motion. "We denote the time-rate 
of variation of h during the motion by dhjdt. From the definition of 
co we obtain 

dh dA dh . 

JTbi +<a bx (120) 

Multiplying by h and using eq 90, we obtain 

Introducing the value of w from eq 112, we obtain 

,dh_ /-rrf Wdh H 2 dhb 2 h hH*& If] 
^dt" i^yjBSat 1 dxdz 2+ 6 dx 3 ]' 

which may be written, after dividing through by h, 

This expression gives the rate of variation of the height of an intumes- 
cence element. 

5. MOTION OF THE CENTER OF GRAVITY OF A WAVE 

The last section treated of the motion of an element, do-, of a wave. 
We may also discuss the motion of a volume, V, between two abscissas, 
Xi and #2 where x 2 ^>Xi. Then 

F- f" do. 

Let £ denote the horizontal coordinate of the center of gravity of 
the volume and rj the elevation of the center of gravity above the level 
of the undisturbed liquid in the channel. By definition, 

rx% rx 2 

I x do h do 

|=^y— ,andi,= ^__. (123 ) 

Differentiating with respect to t, we obtain 

("a da C'-d« 

dt V > &n<i dt 2V ' (1ZV 

which are expressions for the velocity of the center of gravity of the 
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C Xl C x% dh 

volume. The two integrals I 0) dor and I -rrda are readily evaluated 

by means of eq 112 and 122. Writing hdx for d<r, we obtain 

;>w^+if>+f'(§):} 



and 



f^-vK^+WMi)')]:- 



When these expressions are substituted in eq 123 and 124, we obtain 
for the velocity components of the center of gravity 



and 

drj 



-^[S+fiC-(l)')I- o*> 



dt~ 

Equation 125 simplifies if the slopes dA/dx vanish at the end points 
of the intumescence considered, or if the volume V is large, yielding 

(§) 2 =Sr(ff+3„). (127) 

Equation 126 likewise simplifies if applied to the entire body of an 
intumescence. If the intumescence is unlimited, V is infinitely large, 
and the right-hand member vanishes. It also vanishes if the in- 
tumescence is limited, for then both h and dh/dx vanish at the anterior 
and posterior encls of the intumescence. Therefore, the center of 
gravity of an intumescence is maintained at a constant height above 
the primitive free surface, and the velocity of propagation of this 
center of gravity is constant. [5, p. 84]. 

6. VARIATION IN ENERGY OF AN INTUMESCENCE 

It may be inferred from the invariability of 77 that the integral 
fh da, when applied to an entire intumescence, is independent of 
time; that is, the integral 



p 



h 2 dx 



is independent of time. When this integral is multiplied by p#/2, the 
product represents the potential energy of the intumescence, eq 62. 
On the other hand, the kinetic energy of the entire liquid is given by 
the double integral 

fco rH+h 

p/2 dx (u 2 +w 2 )dz. 

JXQ JO 

The velocity components, u and w, are given by eq 114 and 115. 
From these equations, neglecting quantities of small order, we obtain 
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which we multiply by pdx/2 and integrate between the limits x=Xo 
and x— oo. The integral hbhj'dx vanishes at infinity and also, when 
the intumescence is limited, at x . It also vanishes at x when the 
intumescence is unlimited, because this integral is negligible in com- 
parison with the first two terms owing to the length of the intumes- 
cence. Thus the kinetic energy of the wave system is 

The total energy of wave, E, consequently, is 

pgE^pg f " h 2 dx+^ f" h z dx. (128) 

It consists of two terms, of which the first is independent of time, 
and the second which is small in comparison with the other, may be 
shown to vary very slowly with time. 

7. THE MOMENT OF INSTABILITY OF AN INTUMESCENCE 

It was seen that every intumescence is characterized by two 
integrals, E=fhd<r and V—fda^ which are independent of time. 
(See eq 65, 69, and 88). In addition to these, there also exists another 
integral 

which is independent of time. The invariant character of this 
integral was discovered by Boussinesq, who referred to it as the 
moment oj instability [5, p. 97], 
We consider the integral 

M *~w E =h Lw ~ rvJh (130) 

where E is the integral representing the energy of the intumescence. 
Differentiating both members of the equation, we obtain, since 
dE/dt=0, 

dM s _ n/dAV 3A 8 _6A 2 j , bh Nf\ 

dt~- IKbx) ip ST^di'd*!. (131) 



J» bt\ 12 H 2+ EP +2 crff*' 



The term in the brackets vanishes when the intumescence is limited; 
in the contrary case it is negligible in comparison with the integral, 
owing to the great length of the intumescence. Therefore, 



dM, 

dt 



! J*, Dt\ 12 H* + H i+2 dx*/ 1 *' VW 

or, in view of eq 112, 

dM 3 12 C'bh, , , . 
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Now, introducing bhjdt from eq 90, and then integrating, we obtain 

If the intumescence is limited, the right-hand member vanishes: if it is 
unlimited, this member is a quantity of the order which is neglected 
in this analysis. Therefore, to this approximation, the moment of 
instability of an intumescence is independent of time. 

8. SOLITARY WAVES 

The effect of height and curvature on the deformation of waves 
was discussed in the previous section on the basis of the assumption 
that deformation actually occurs. It may well be that configurations 
exist with appreciable height and curvature which are propagated 
with no deformation whatever. As a matter of fact, such waves were 
first discovered by experiment. The solitary wave of J. Scott Russell, 
consisting of a single elevation whose height may be of the same order 
of magnitude as the depth, may travel considerable distances without 
appreciable deformation [6, p. 505; 7, p. 324; 9, p. 289]. This property 
is called longevity and the law governing it may be derived from eq 106. 

Putting new variables 

(a:— o#)=rH, 

h= v H, 

and 



G4-0- 



b 



b 2 d 2 

into eq 106, we get, since 7^2 ==aj2 3r2 on the basis of the assumption 

that the wave does not deform, 

^ 6677 _ V _ 2 p^ 0j (135) 

which when integrated twice gives 

667 7 -9) ? 2 -2^+c 1 f+c 2 =0. (136) 

Since rj and ^-^ are finite, ^=0. If we restrict ourselves to the 

52 
case where at £= 00 f 77=0 and =^72=0, then c 2 =0 and we have 

6677— 9t7 2 — 2^=0. (137) 

dw 
Multiplying by ~ d£ and integrating: 
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(138) 



Denoting the maximum height of the wave by hi, and the corres- 
ponding value of 7} by rj ly 






(139) 



At the crest of the wave, ij=r) U and the wave surface is parallel to the 
bed, d?j/d|"=0. Substituting these values in eq 138 we obtain the 
relation between b and ijj, 



(140) 



(141) 



Replacing this value of b in eq 138: 

(!)'=3„'„-3,., 

and integrating twice: 

^sech{^-} 

Reverting to the original variables, we obtain the equation of the wave 
profile 



h=h t sech 2 ! J ||p(3— «*)!• 



(142) 



The free surface is symmetrical about the normal at x=cat, and 
approaches the primitive surface as x— cat approaches ± oo (fig. 2). 




*• 4- x 



Figure 2. — Diagram of a solitary wave. 



The velocity of propagation, <o, which may be interpreted either as 
that of the volume element or height element since the wave does 
not deform, is obtained from eq 140: 



J-giH+h), 



(143) 



which is the formula adopted by Russell [7, p, 328; 8, p. 423], on 
empirical grounds and confirmed by the experiments of Bazin [6, p. 
496, 515]. 

198881—39 6 
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A comparison of theory and experiment, based upon 117 observa- 
tions by Bazin [6, p. 510-514], shows the error of the formula to be 
— 0.91 ±5.10 percent. Similar calculations based upon 77 observa- 
tions by Russell [8, p. 440-442] give +0.78 ±2.30 percent. The ± 
term indicates the standard deviation, and this usage will be followed 
in subsequent comparisons. It will be noted that in each case zero 
error is well within the limits determined by the standard deviation. 

From eqs 127 and 143 we will now show that the center of gravity 
of a wave propagated without change of form is Z^/3 above the primi- 
tive surface, where hi is the maximum height of the wave. Since the 
wave retains its form, 

|= W (144) 

and is constant. Substituting from eq 143 in eq 127, we obtain 

U=|- (145) 

A relation connecting the height at any arbitrary point with the 
partial volumes of the wave on each side of the perpendicular at this 
point is derived as follows. By virtue of the relation(l/A) (d/dx) = d/dc, 
when dhjdt is placed equal to zero, eq 122 yields the result: 



„/, . 2H*d 2 h\ . . 

hH 1 H — o~- xj )= constant. 



The constant is zero, since h vanishes at infinity; therefore, 

,.2H*d 2 h n nAa . 

1+ T" ^ =0 - (146) 

Integrating and employing the boundary conditions <r=0, h=Q; 
<r—Q, ft— 0, we obtain the solution 

h=^{Q-e)e. (147) 

Since the wave is symmetrical, the maximum ordinate, h i7 divides 

the wave into equal parts: r)=h 2 , c=Q/2. 

Therefore, 

Ai = jHb> (148) 

which is the relation between the maximum height of the wave and 
its total volume. It was demonstrated in section IV-2 that the 
potential and kinetic energies of an intumescence are equal, to within 
a small term of the order of h z ; consequently 

E=2r,Q=~^. (149) 
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From eq 148 and 149 we derive the following relations: 

*»— &T faad 2-8^ff (150) 

These relations will now be used to explain the phenomena of waves 
traveling up a channel of slowly decreasing depth. Considering eq 
85, we see that the quantity Ci—c 2 is small, since dll/dx is small. 
Consequently, the reflected wave is negligible, and the energy of the 
transmitted wave remains sensibly constant. From eq 150 we deduce 
that hi continually increases, and Q/Ih continually decreases. The 
wave gets higher and higher, and its base narrows continually. As 
a result, the wave becomes unstable and finally breaks. A quantita- 
tive analysis of this phenomenon will be given in section V-19. 

9. CNOIDAL WAVES 

In the theory developed thus far, extensive use has been made of 
the assumption that the liquid is undisturbed at infinity ; that is, that 
the wave profile x=j(h) vanishes at infinity together with all its 
derivatives. In the previous section this assumption was utilized in 
deriving eq 137, and it was found that the only wave fulfilling this 
condition and propagated without change of form was the solitary 
wave. We now wish to eliminate this restriction and determine the 
form of a wave of appreciable height, governed by eq 106 and propa- 
gated without change of form. 

We therefore assume that atA=0, the derivatives ph/dx, b 2 h/bx 2 f 
... do not all necessarily vanish. The constant c 2 is consequently 
retained in eq 136, and we write 

66if-V-p+3c2=0. (151) 

Multiplying both members of eq 151 by (drj/dt;) df and integrating: 

(^=Mrf-W+Zc 2 r)+Zc u (152) 

where 3c i is another constant of integration. Since drj/d^ is real at 
77=0, then C\ is positive, and if rj' ', 77", 77"' are the roots of the poly- 
nomial, then 

V V V — Ci. 

Either all the roots are positive, or one is positive and two negative, 
or one is positive and two complex. There is at least one positive value 
of 77, say 77', for which drj/d'f vanishes. If this is the only value of 77 for 
which the profile is horizontal, then the wave profile must become 
infinite with f , which is physically impossible. The third case is there- 
fore excluded. If all the roots are positive, then the wave profile must 
be included between the extreme values, since otherwise it would 
become infinite. This requires that the wave surface be everywhere 
above the primitive surface, which is likewise physically impossible, 
and the first case is also excluded. Therefore 77" and 77"' are both 
negative. Writing the three roots ^,—772,-773, 773^*72, we have 



~~ 



82 



and 
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<WJ 



= 3 (vi — v) (v + V2) (v + fa) , 



C2 = — Wz + ViVz + ?7i*?2 j 
Ci — 7jiTj 2 ris. 

Introducing the new variable, x, 



eq 153 becomes 



where 



rj=r jl cos 2 X—V2 sin 2 x , 



/?(^)=Vl-^ 2 sin 2 x, 



and 



P \nrii+Vz)~~\Hhi + h 



+hr 



_Vi+V2_h 1 +h 2 



(153) 
(154) 

(155) 
(156) 

(157) 
(158) 



Vi+Vz h+h 
Selecting f =0 at the apex, then y=r}i and x=0, and we have 

r-* r /1 f 2 x • 2 = ^(*> *), (159) 

Jo Vl— * sm x 

where F is an incomplete elliptic integral of the first kind. Introduc- 
ing Jacobian elliptic functions, since cos x=cn(f//3, k), sin x— 
sn(f//3, k) and sn 2 w+cn 2 u=l, we have 



-(*-«*)' vd^j 



(160) 



(161) 



or reverting to the original coordinates 

This equation represents an infinite number of undulations of 
identical size and shape, each symmetrical about a vertical plane 
passing through the apex. The existence of these periodic waves 
moving without change of form was first indicated by Boussinesq 
[10, p. 392], and the name, cnoidal waves, from the en function by 
analogy with sinusoidal, is due to Korteweg and de Vries [11]. 

The relative wave length, A/fiT, is given by 

^=20 P /2 & , , =2pF l (k)f (162) 

ti Jo VI— Arsnrx 

Fi being the complete elliptic integral of the first kind. 
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The initial depth of the water in the channel being H, we must 

have 

x 



=0, 

or from eq 155 and eq 156 



J o H wr=( 



i 



2 fa cog 2 x—V2 sin 2 x)J d ==Qj 



-y/l—k 2 sin 2 x 

which from eq 158 reduces to 

(h l +h a )E l (k)=hjr 1 (k) 9 (163) 

where Ei is the complete elliptic integral of the second kind. The 
form of the waves is completely determined by \/H, and hi/H, the 
relative wave length and height of the crest above the initial surface 
of the water, since h 2 /H and hJH are determined by eq 158, 162, and 
163. 

From eq 162 we see that the wave length of the undulations be- 
comes infinite when k=l. From eq 163 if fc=l, h 2 —h z =0. Since 
cn(x, l) = sech x, then for &=1, eq 161 reduces to the equation of the 
solitary wave (eq 142), and the latter appears as a special case of the 
cnoidal wave. 

From the definition of b (eq 136) and from eq 154, we have 

j^+BTTTT (164) 

which is the expression for the velocity of the wave profile. Since 
we are now dealing with an infinite number of undulations moving in 
a channel of finite depth, a question arises as to the meaning of the 
velocity of propagation. Stokes [12, p. 202], has given two definitions 
applicable to this case, and we adopt the second one; namely, the 
velocity of propagation is the velocity of the wave profile with respect 
to the center of inertia of the fluid mass included between two vertical 
planes separated by a distance equal to the wave length. 
Writing this velocity as a> c : 

co c =co-^ J U(H+h)dz, (165) 

where U is the mean velocity in a vertical section. If U is the mean 
velocity in the section at h= 0, we have from eq 95 

U(H+h)=hcc+ UoH. (166) 

Substituting in eq 165, we have for the velocity of propagation 

o, c =co-?7o (167) 



or 



fcV 



! =* 2 (l-2-^) (168) 



We now evaluate U Q . Integrating eq 114, we have 

TJ H 2 &u . n 
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u being the velocity at the bottom. It has been assumed in the 
analysis that at h=Q, the surface velocity u x vanishes (see the discus- 
sion of eq 41); hence from eq 114, 



and 






(169) 



From eq 104, since we can write b/Z)t= -co (d/bz), we obtain by differ- 
entiating twice with respect to x, and neglecting the term containing 
d%/dx 4 : 

Substituting in eq 169, we finally have 



and eq 168 becomes 






or from eq 151, 154, 164, in terms of h u h 2 , and h z , within the approxi- 
mation of the analysis, that is, neglecting h^Th/SPy 



w;= ,h[i + |][i-|][i-|]. 



(170) 
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Figure 3. — Cnoidal waves. 
Relative displacemen of the free surface as a function of the fractional part of the relative wave length. 

This, then, is the velocity of propagation of the wave profile with 
respect to the center of inertia of the fluid mass included between 
two vertical planes separated by a distance equal to the wave length, X. 
For other treatments of this subject reference may be made to 
Boussinesq [10, p. 390], Korteweg and de Vries [11], and Lamb [4/p. 
402]. 
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Figures 3 and 4 represent two typical cnoidal waves. Figure 3 
shows on an arbitrary scale one complete undulation for two values 
of X/jEZ, and figure 4 shows extended profiles of the same two waves 
at the correct height above the bottom, but compressed horizontally 
5 times. In both cases, hi/H—0.25. 
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Figure 4. — Cnoidal waves. 
Relative height above the bottom as a function of relative horizontal dimension. 

10. MOMENT OF INSTABILITY OF A SOLITARY WAVE 

Of all the intumescences of equal energy, the solitary wave is the 
only one for which the moment of instability 



M ^Ll(W^ dx 



(129) 



is a minimum. 

The theorem may be proved by the methods of the calculus of 
variations for the solution of isometric problems. By an isometric 
problem, we mean one of the following kind: 

Among all curves joining two points (x 0} h ) and (x h hi) for which 
the definite integral 



K 



= P 1 G(z, K x', h')dx 



takes on a given value I, to determine the one 

Co, x=x(s), Xq=x(s )\ 
h=h(s), ho=h(s )\ 

which minimizes (or maximizes) another definite integral 



(171) 



(172) 



J= 



Fix, h, x\ h')ds. 



(173) 
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Here s is a parametric variable, and the primes indicate differentia- 
tion with respect to s. 
Putting 

H=F+\G, (174) 

where X is a numerical constant, the first necessary condition for a 
minimum or maximum is that the curve C in eq 172 must satisfy 
the differential equation 

H xh ,-H x , h +H 1 (x'h"-x"h')=0, (175) 

where 

The subscripts indicate partial differentiation with respect to the 
variables represented by the subscripts. The second necessary 
condition for a minimum is that 

#i^0 (177) 

at every point of the curve C . 

Obviously, in the present isometric problem, the integral K corre- 
sponds to the energy of an intumescence and the integral J to the 
moment of instability; that is, 

K= f °V<fo, (178) 

J xt 

and 

^r[(S)-i>. 

We select the parametric variable s so that 

x=s, ] (180) 

and } 

h=h(s).\ 

The above integrals, eq 178 and 179, become 

K= (°h?x'ds, (181) 

and 

J=r(^—ah?x'\ls, (182) 

where a=3/IP. Hence, 

G=¥x', (183) 

F=^—ah 3 x', (184) 

fli-J* (185) 



and 

H -x 



H=H-cth?x'+Wx'. (186) 
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The first necessary condition for a minimum, eq 175, after some 
simplifications, yields 

Zah 2 -2\h+2^=0, (187) 

which also may be written as 

3 . Iff VA_1 ™ ( . 

if a is replaced by 3/H 3 . Since the left-hand member is a constant, we 
see on comparing eq 188 with eq 112 that the intumescence having a 
maximum or a minimum moment of instability for a given energy is 
one that is propagated without change of form. The only wave which 
is propagated with constant velocity is the solitary wave. The second 
necessary condition, eq 177, is likewise satisfied; hence the moment 
of instability is a minimum. 

11. MOMENT OF INSTABILITY AND WAVE FORMATION 

The frequent occurrence of the solitary wave, owing to its produc- 
tion by a broad class of different initial disturbances, was observed 
by Scott Russell [7, p. 351; 8, p. 423]. This frequent formation of the 
solitary wave has been explained by Boussinesq [5, p. 99; 10, p. 401], 
utilizing the theorem just proved — that the solitary wave is the only 
intumescence for which the moment of instability is a minimum. An 
integration would show that this minimum is 

97 77W 

(M.)i=$£g[, (189) 

where E is the energy of the wave and H is the depth of the channel. 

If one considers an intumescence of arbitrary form, and a solitary 
wave, both having the same energy, E, the excess of the integral M s 
of the intumescence over that of the solitary wave, M s -~ (M t )i may 
be regarded as a measure of the rapidity with which the intumescence 
will deform during its propagation and also of the extent of the 
deformation taking place. This is why Boussinesq named the integral 
the moment of instability. 

If M s — (M s )i is small, the profile of the intumescence will oscillate 
about that of a solitary wave having the same energy, E. The form 
of the intumescence cannot change very markedly, as this would 
require M t to increase, contrary to the theory. Friction, which has 
been ignored in the theory, will damp out these oscillations and the 
intumescence will eventually assume the shape of a solitary wave. 

It should be noted that there does not exist a solitary wave having 
an arbitrary energy, since there is a maximum height, above which 
a solitary wave cannot be propagated without breaking (see section 
V-19). This limiting factor serves to determine the class of disturb- 
ances which can produce a solitary wave. 
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12. NEGATIVE WAVES 

Let us now examine the circumstances governing the motion of a 
negative wave limited in extent and propagated in the positive 
^-direction, negative signifying that the surface is depressed. The 
longitudinal profile of such a wave in its simplest form will consist 
at the moment of the formation of the wave of a main concave por- 
tion in the center and two convex arcs at the posterior and the anterior 
ends of the wave. As long as the form is not altered, the center of 
gravity of the wave is located at one-third of its maximum depth; 
that is, rj=hif'3. According to eq 127, the velocity of propagation of 
this wave would be 

J^giH-ht), (190) 

a law that is confirmed very nearly in the experiments of Bazin. 
Calculations based upon 18 observations, [6, p. 521] give as the error 
of the formula, —1.54 ±2.09 percent, where the ± term indicates 
the standard deviation. 

The initial simple form of the wave, however, is soon modified. 
In fact, eq 112 shows that since h is negative, the third term is nega- 



//Xnx 
Figure 5. — Negative wave as observed by Bazin, 

tive in the concave portion and positive in the convex portion. The 
second term is more negative in the concave portion than in the con- 
vex portion, since h is smaller in absolute magnitude in the latter 
case. Thus the head of the depression, being both at a higher eleva- 
tion than the center of the cavity and convex, while the latter is 
concave, is propagated faster than the center and becomes elongated. 
The tail of the wave also travels faster than the center for the same 
reasons and is therefore shortened as long as the requisite initial con- 
ditions obtain. However, the process of shortening cannot be indefi- 
nite. Furthermore, a convex portion is necessary for the continuity 
of the surface slope ; hence the convex portion is finally raised above 
the level of the primitive free surface. Since h is now positive, the 
last term of eq 112 is negative; and if either h is small enough or 
d 2 hj dr* large enough, it can exceed the second term, so that a condi- 
tion is reached such that w of the elevated convex portion is the same 
as co of the lowest point of the depressed body. But the elevation 
thus produced at the posterior end of the main negative body of the 
wave requires an intermediate concave portion in order to join the 
free primitive surface with a continuous slope. Such a concave por- 
tion cannot exist and cannot be formed if its velocity of propagation 
is greater than that of the positive wave preceding it. Equation 112 
indicates that if these velocities are to be equal, the convex positive 
wave must be followed by a negative wave. Continuing the argument, 
we see that a given negative wave is followed by a series of smaller 
waves, alternatively positive and negative, and entirely convex and 
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concave, with inflection points on the level of the primitive surface 
(fig. 5). For additional figures, see [6, plate 2, fig. 3; 7, plate 52, 
fig. 10]. 

13. INTUMESCENCES OF FINITE HEIGHT BUT NEGLIGIBLE 

CURVATURE 

For an intumescence of finite height but having a free surface of 
negligible curvature, the expression of the velocity of propagation 
of volume element, eq 112, assumes the simple form 



= V<7#(l+f£)- (191) 



As an immediate consequence of this simplification, differential 
equation 90 can be integrated and the equation of the intumescence 
profile obtained. Introducing eq 191, eq 90 becomes 



&H i+ lm=°- <i92) 

which possesses the integral 

x-jm(l+%%y=M), (193) 

where jih) is an arbitrary function of h. The integral may be inter- 
preted as follows: If h is the elevation of the free surface at x at time 
t y then the same elevation is reproduced at x+Ax at time t+At; but, 
from eq 193: 



A*=VffS(l+§^)«, 



and, therefore, the intumescence height is propagated with the 
velocity 



«i= 



s Vfffl(l+§g). (194) 



which is Airy's expression [4, p. 243] for the velocity of propagation 
for waves of finite height. It should be noted that eq 194 represents 
the velocity of propagation of the height of an element. 

The arbitrary function,/, in eq 193 can be determined in two ways, 
depending on whether the canal is infinite in one or both directions. 
In the latter case, h is usually given for all values of x at some instant, 
t. If the relationship is in the form h=fi(x), we may invert the 
function and obtain f 2 (h)= x, which is the desired function. 

In the former case, assume the channel in communication with an 
infinitely large basin and that the water level of the basin is a known 
function of time, say, h—F(t). Inverting this function, we have 
t=Fi(h) } and, if we take the point where the channel meets the basin 
as the origin of x } the arbitrary function is 



/(A)=-V^(l+|^i(A). 
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In intumescences of this type, the mean velocity, U, of the fluid 
particles in a cross section is 

17-VS». J(l-Ej) (195) 

which follows from eq 96 and 191. 

14. HORIZONTAL CHANNEL IN COMMUNICATION WITH THE OCEAN 

Boussinesq has applied the method of the last article to a horizontal 
channel in communication with the ocean to determine the mean 
elevation of the ocean surface above the bed of the channel. 

We may write for the elevation of the surface water of the ocean 
above the bed of channel, 

H+h=a(l + a sin ?~ \ s=0, (196) 

where H is the depth of water in the channel at points so remote from 
the mouth that the disturbances of the ocean are not felt, a is the 
elevation of the mean level of the surface waters of the ocean above 
the channel bed, a the ratio of the half amplitude of the tidal oscilla- 
tions to the elevation, a, and T the period of the tides. Using eq 
171, we obtain for the arbitrary function, f(h), in eq 193, 

-/«-v*a(i+§&)£*r* f^> (197) 

and for the free surface of water in the channel, 

-vs<>+iiX-£-*' £ S= s > w 

After a certain length of time, the transients will die out and the 
total volume of water entering the channel during each period T will 
be zero; that is, 

f T U(H+h)dt=0, x=0. (199) 

This condition permits the determination of H. Since U(H+h)—o) h, 
eq 199 gives, substituting for o> from eq 191, 

Introducing h from eq 196 and integrating, we obtain 

„ 3 {H-aY _Zc?a 2 

M a 4 H ~8 H ' 

which can be simplified further since H—a is a small quantity: 

ff-a=! a 2 a. (201) 






ftitoSft] Irrotational Translation Waves 91 

This is the desired relation. The mean level of the surface water 
of the ocean is below the level of the water in the horizontal channel 
at large distances from the mouth by the small amount (3/8) a 2 a. 

15. WAVES DUE TO SUDDEN INCREASE OF DISCHARGE 

Let us examine the case in which a positive intumescence of infinite 
length is produced in a body of still water in a channel when a dis- 
charge of constant value is established at the entrance. After the 
initial disturbances have died down, the wave is of constant elevation, 
h, except in the anterior portion, and the front of the wave moves 
with a constant velocity, co. As the curvature of the layer represent- 
ing the intumescence is negligible at points sufficiently removed 
from the front, the velocity of propagation of the wave front must 
be from eq 191, 

»=y ff(flr+| a), (202) 

since h/H is small. 

Bazin's experiments verify this formula closely. From his data 
[6, p. 539-541], velocities were computed and were compared with the 
theoretical values. The mean values of h and H at the two points of 
the channel between which the wave traveled when its elapsed time 
was observed were substituted in the formula. The method of 
obtaining h from the data is described on page 96. 

The data gave sufficient information to make 37 comparisons. 
Observations on eight intervals were rejected, since the observers 
reported that the wave had completely broken and reported h as zero. 

The average error of the theoretical values was — 0.98±2.71 per- 
cent, where the ± term indicates the standard deviation. These 
values were plotted as a function of h/H in figure 6, in an attempt to 
discover a systematic variation. The lines connect observations on 
the same wave in different intervals. The results were essentially 
negative, although there is a suggestion of an error changing in sign 
from minus to plus, as h/H increases. It will be noted that the for- 
mula applies equally well to unstable waves on the point of breaking 
as well as to those having a definite hijh ratio of 3/2. (See fig. 9.) 

St. Venant, using the principle of momentum, has derived for the 
celerity of this type of wave, 

W 2 =^(l+l^+i|) (203) 

an expression which reduces to eq 202 when the last term is neglected. 

16. EFFECT OF CHANNEL VELOCITY DISTRIBUTION 

We may now consider a discharge wave superposed on a current 
and determine the effect of the velocity distribution in the current on 
the celerity of the wave. 

We shall assume (1) that the profile of wave does not deform dur- 
ing the motion and (2) that at the two sections, each at a distance s 
from the wave front, the velocity distributions are characteristic of 
the channel. (See fig. 7.) Let the mean velocities at these sections 
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Figure 6. — Percentage of error of eq 202 compared with experiment, as a function 

of the relative height. 
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be Ui and U 2 , the local velocities at any height z be u x and u 2 , and the 
depths H and H+h. We consider the body of water included at the 
instant t between the sections AB and Du. The distance between 
AB and the section containing the wave front is s; the distance be- 
tween DC and the wave front is 8+coAt, where a is the celerity of the 
wave, and At an interval of time. The momentum of the body of 
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Figure 7. — Diagram illustrating the effect of channel velocity distribution. 
(a) Configuration at time til (&) Configuration at time h+At. 



water thus delimited may be denoted by M x +M 2 , where M 2 is given 
by 



M 2 =pwAt I u 2 dz=pa)AtU 2 H. 



At the instant of time t-\-At, the particles originally at AB have 
reached the position A'B\ and the particles originally at DC have 
reached the position D'C' . (See fig. 7.) The liquid contained 
between A'B' and D r C f has the momentum M x +M\+M' 2 , where 






(uAt— UiAijUi dz, 



and 



M' 2 =pj o 



u 2 At dz. 



Introducing the Boussinesq number , 77, the coefficient of the excess 
of the mean of the velocities squared over the square of the mean 
velocity, which we assume to be unaffected by depth, 
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X 



ufdg^a+vWiH+h), 



and 

("u t 'dz=(l+r,)U a 2 H. 

We have, therefore, 

M;= P Z7i(ff+A)«A*— pU, 2 (H+h) (l+ri)At, 
and 

M,=p(l+r,)U a * HAt. 

Thus the increase of momentum of the body of water under con- 
sideration during the time interval At is 

Mt+MS+M,'- (M x +M 2 ) = P At[U x (H-\-h)(co- U t ) + 

U 2 H(U 2 -co)- v { U x 2 {H+h)-U 2 2 H}]. (204) 

We may eliminate U\ by the equation of continuity, 

U 1 (H+h) = U 2 H+o>h, (205) 

and obtain for the rate of change of momentum: 

H+h l 1+ i 1 —^-U^]j> (206) 

neglecting h/H with, respect to unity. 

The force which produces this rate of change of momentum is the 
difference between the total hydrostatic pressures acting on the 
sections AB and DC, which is 

P0fl&[l+gg]- (207) 

Equating 206 and 207, and assuming that the effect of the 77 term 
is small, we obtain 

or 

-^^(i+i^x^C 1 -?^ 5 ]) (208) 

This is the expression for the celerity of a discharge wave or a de- 
pression wave, which takes into account the height, h, and u 2 , the 
channel velocity distribution. For a wave moving downstream the 
plus sign holds; for a wave moving upstream the negative sign holds. 
For a positive wave h represents an elevation and is positive; for a 
negative wave h represents a depression and is negative. 

Equation 208 can be further simplified when the approximation 
w=Z7 2 ±Vfl f S' i s used. Accordingly, we obtain for a wave moving 
downstream 
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„= K -V^(l+|^)(l+l[||+^]). W 

and for a wave moving upstream 

«=^-v,«(.+^)(. + f[i-^])- *m 

It is now seen that the effect of velocity distribution is appreciable 
only when the ratio UflgHia larger than unity, and for the ordinary 
case, in which U 2 2 /gH is a small fraction, this effect is negligible. 

17. INITIAL WAVES 

Bazin observed that in waves of constant discharge moving in still 
water, the initial part of the wave is made up of a series of relatively 
short, alternately convex and concave, intumescences. (See fig. 8.) 



Figure 8. — Initial wave as observed by Bazin. 

The foremost intumescence is convex, and its crest is at a higher lever 
than that of the constant elevation of the free surface of the main 
intumescence. Bazin refers to this initial part as the initial wave. 
The mechanical basis of the phenomenon was discussed by Boussinesq 
in the following manner, which is similar to the treatment in section 
V-12. 

At the head of the intumescence it is not possible for the depth to 
equal the depth of the main layer. The anterior part of the wave 
being convex, the curvature i> 2 h/dx 2 is negative and according to eq 
112, 

"=^\ 1+ lH+Mw) f (112) 

the velocity of propagation there is less than that part of the layer 
which follows. Tins layer moves with a greater velocity and thus 
overruns the anterior portion of the wave until the augmented height 
compensates for the effect of the negative curvature. This explains 
the formation of the initial wave. It is known that the height of the 
initial wave tends to oscillate around a height equal to 3/2 of the 
height, A, of the long layer. The height of the initial wave will increase 
until the velocity of propagation of its center of gravity becomes 
sensibly equal to that of the layer. But the form of this wave indi- 
cates that the center of gravity, rj, is to be placed at 1/3 of its height, 
h a ; and, on the other hand, eq 127, which is applicable in the present 

198881—89 7 
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case, gives for the square of the velocity of the propagation, g(H+h a ). 
Equating this to the square of the velocity of the layer, g(H+3hj2), 
eq 202, it is seen that h a =3h/2. 

The ratio hi/h=3/2 has been investigated by Bazin [6]. The ex- 
periments were made in a rectangular channel having a bottom slope 
of 1.5 in 1,000. The percentage change in depth being the least at 
the deep end, where the waves were generated, the motion at this end 
will nearly approximate flow in a, channel of uniform depth. Calcu- 
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-Relative height of the head of an initial wave as a function of the relative 
depth of the main layer. 

lations made by the authors from data taken at his end of the channel 
[6, p. 542] give Ai/A== 1.61 ±0.18. Bazin [6, p. 550] considers 3/2 
to be the approximate value. 

Inasmuch as the waves are found to break as they ascend the 
channel if the original depth of the water is small enough, it is instruc- 
tive to consider this aspect of the phenomenon. Observations were 
taken on the wave at four points in the channel. The complete 
data [6, p. 539-541] have been employed to plot hjh as a function of 
the dimensionless ratio h/H, the relative height of the wave layer. 
The average value of h was utilized in order to approximate the 
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steady height of the layer, eliminating the undulations behind the 
head of the wave. These data are represented in figure 9. The 
circles represent observations on waves which have no tendency to 
break; the squares, those in which breaking is imminent; the triangles, 
those in which breaking has occurred. An average of the 33 points 
with circles gives Ai/ft=1.51±0.17. Although the data do not pro- 
vide a clearly defined function, the following conclusions seem justified. 

1. As would be expected, for increasing h/II the waves become 
unstable and eventually break. The critical value is about 0.275. 

2. The value Ai/^=1.5 is a reasonable approximation in the stable 
region. 

3. The ratio hi/h definitely decreases with increasing h/H, outside 
the stable region. 

Recent experiments [2, p. 158-159] indicate that hi/h is a function 
of h/H in the stable region, increasing from 1 to 2 at the critical 
value 0.28, then decreasing asymptotically toward unity in the 
unstable region. 

Where the initial wave joins the main layer, the condition of con- 
tinuity of slope requires a concave curve, that is, a curve for which 
d 2 h/dx 2 is positive. The positive curvature would entail a velocity 
of propagation greater than that of the layer as a whole, unless h is 
less than in the remainder of the layer. Hence a depression is pro- 
duced immediately following the initial wave. 

Continuing the argument, it is seen that the initial wave is followed 
by convexities and concavities all above the undisturbed surface. 
Because of the internal friction of the liquid, the heights of these 
convexities and concavities decrease in the posterior direction until 
they become insensible. The inflexion points where d 2 h/dx? vanish 
are in a plane containing the free surface of the main layer. 

18. FAVRE'S EXPERIMENTS ON THE UNDULATIONS IN A POSITIVE 

SURGE 

In 1934, Favre [2] made an elaborate series of experiments on trans- 
lation waves. Several of these experiments were undertaken to 
investigate the undulations composing the head of a discharge wave 
[2, ch. IV]. The experiments established the following facts: 

1. The undulations are not formed immediately, but require a 
certain amount of time to establish themselves. 

2. After the undulations are formed, they are of similar size and 
shape. 

3. The final configuration is a stable one. 

4. The height of the undulation at x is independent of 9, the time 
in which the flow is established, for a considerable range of values of 
x, provided 9co/x<0.2. In other words, within certain limits, a dis- 
charge which is not created instantaneously will behave like a theo- 
retically sudden increase of flow. 

5. The wave length of the undulations depends only on the mean 
height of the layer above the primitive surface. The former decreases 
as the latter increases. 

These results suggest that the undulations composing the head of 
a discharge wave have a very general character, in the same sense 
that the solitary wave has a general character. Within certain limits, 
they are formed independently of the rate at which the discharge is 
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increased, and depend only on the final discharge rate; they assume 
a stable configuration after a short initial interval, that is, they are 
propagated without change of form. In section V-9, it was mathe- 
matically demonstrated that the only wave propagated without change 
of form is the cnoidal wave, which includes the solitary wave as a 
special limiting case. It is also true that the solitary wave is a special 
case of the discharge wave from the physical point of view. The 
solitary wave may be formed by a sudden increase of discharge which 
soon decreases to zero, while the discharge wave is caused by a sudden 
increase of discharge which is maintained. All these considerations 
inevitably suggest the question: Are the undulations composing the 
head of a discharge wave cnoidal waves? 

A comparison of figure 4 with Favre's figure 48 [2, p. 156] supports 
an affirmative answer. The wave profiles are very similar. However, 
there is a more significant comparison which may be made, namely, 
one between the theoretical velocity of propagation, « c , and that 
observed by experiment. This is accomplished by using eq 167 and 
170. From the values of n x and \/H given by the experimental data, 
we determine r; 3 , k, and y 2 by eq 163, 162, and 158, and then calculate 
«,. from eq 170. Equation 191 written in the form 



-V^(i+fzp) 



(191) 



is used to compute «. This equation, in the form of eq 202, has already 
been verified in section V-15. In order to calculate Uo, we write 



wh'=U (H'+h'), 



(211) 



where h f is the average height of the layer, and IT is the undisturbed 
depth of the liquid. It should be noted that this H' is different from 
the II in eq 170. U is computed from eq 211 with the assistance of 
eq 191. From eq 167 we have o) c =a)— U . These two values are 
compared in the following tables for nine nonbreaking waves. Table 1 
gives the experimental data from Favre [2, p. 157], and table 2 gives 
the results of the calculations. The explanation of the symbols is 
clearly indicated in figure 10. 




h' Uo=0 



Figure 10. — Diagram illustrating the notation used in the studies of the positive 

surge. 
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Undis- 
turbed 
depth, 
H' 


Mean 
height 
of layer 


Crest height, h*, 


Mean 
h* 


Wave 
length 




Experiment No. 


1st 
wave 


2d 
wave 


3d 
wave 


4th 
wave 


h'!H' 


4 


cm 
20.55 
20.22 
20.19 

20.38 
20.57 
10.78 

10.73 
10.79 
10.74 


cm 
2.29 
3.08 
4.08 

4.84 
5.72 
0.86 

1.50 
2.48 
3.02 


cm 
3.90 
6.00 
8.10 

9.80 
11.70 
1.30 

2.80 
4.95 
6.25 


cm 
3.60 
6.00 
8.00 

8.70 
11.40 
1.20 

2.70 
4.70 
6.10 


cm 


cm 


cm 
3.75 
5.97 
8.05 

9.47 
11.33 
1.14 

2.76 
4.70 
6.04 


cm 
222 
194 
192 

185 
176 
119 

101 
91 

86 


0.111 


«* 


5.90 




.152 


8 


.202 


10 


9.90 
10.90 
1.00 

2.85 
4.70 
5.70 


"T65" 

2.70 
4.45 
6.10 


.238 


12 


.278 


21 


.080 


22 


.140 


23 


.230 


24 _ 


.281 







Table 2. — Velocity comparisons between cnoidal waves and positive surges 



Experi- 
ment No. 


Ai=mean 
h*-h' 


H=H'+h' 


\{H 


m 


W 


VI 


sin- 1 k 


Oi 


U* 


u—Uq 


b>e 


Per- 
cent- 
age of 
differ- 
ence 


4 


cm 
1.46 
2.89 
3.97 

4.63 
5.61 
0.28 

1.26 
2.22 
3.02 


cm 

22.84 
23.30 
24.27 

25.22 
26.29 
11.64 

12.23 
13.27 
13.76 


9.72 
8.33 
7.91 

7.34 
6.69 
10.22 

8.26 
6.86 
6.25 


0.064 
.124 
.164 

.184 
.213 
.024 

.103 

.167 
.219 


0.052 
.094 
.118 

.126 
.157 
.022 

.081 
.129 
.166 


0.142 
.197 
.222 

.257 

.307 
.126 

.198 

.289 
.349 


48.6° 
55. 4° 
58.8° 

58.1° 
57.5° 
33.6° 

51.5° 
53.7° 
55.4° 


cm/sec 
153.8 
156.9 
162.0 

166.7 
171.7 
109.0 

113.4 
120.6 
124.3 


cm/sec 
16.4 
20.7 
27.2 

32.0 
37.4 
8.0 

13.9 
22.5 
27.3 


cmfsec 
138.4 
136.2 
134.8 

134.7 
134.3 
101.0 

99.5 
98.1 
97.0 


cm/9ec 
139.2 
136.7 
137.9 

137.9 
135.2 
100.0 

98.7 
97.0 
94.5 


+0.58 


6 


+.37 


8 


+2. 30 


10 


+2.38 


12„._ 


+0.67 


21 


-.99 


22 


-.80 


23 


-1.12 


24 


-2.58 







The average percentage error of co c as compared to a)—Uo is +0.09 
±1.53, where the ± term indicates the standard deviation. This 
is further evidence that the undulations in the head of a discharge 
wave are cnoidal waves. Inasmuch as these experiments were not 
intended to shed light on this particular problem, other points of 
verification cannot be easily attempted. An experiment made especi- 
ally to investigate this problem would be more suitable. 

19. THEORY OF BREAKERS 

As was indicated at the end of section V-8, there is a limiting 
height for a solitary wave, dependent on the depth of the liquid in 
the channel. If an attempt is made to form a higher wave, it breaks 
at the crest. If co is the velocity of propagation of the solitary wave, 
and Ui is the absolute velocity of the liquid particle at the apex, that 
is, at the point h=h h it is evident that the limiting height is repre- 
sented by the relation 

Ul — o)=0, (212) 

since if ^i>o>, the particle at the crest of the wave is traveling faster 
than the wave itself, and consequently the wave is breaking. 
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Equating eq 112 and 143 and neglecting terms in (hi/H) 2 or higher 
powers, we obtain for the curvature at the apex of the wave 



/b 2 h\ ^^3 



3 /t! 2 



Substituting this value in eq 116, putting z—H+hi and neglecting 
terms of the order of (h JH)* } the resulting velocity of the particle is 

WHh+!!)} (2i3) 

Putting this value and that from eq 143 in eq 212 ? the limiting 
height is a root of the equation 



!+K!)W 1+ b- (2I4) 

The only positive real root of this equation is h^H =0.731 } and this 
is the limiting height of the solitary wave. 

This figure is confirmed reasonably well by experiment. Calcula- 
tions made by the authors from Bazin's data [6, p. 510-514] give 0.71. 
Scott Russell [7, p. 352; 8, p. 425, 426, 445, 450] surmised that the 
wave broke when its height was approximately equal to the depth, 
and McCowan [13, p. 56] states that for his experiments 0.75 is a 
closer approximation for the relative elevation at the breaking point. 

The above result may now be applied to the ultimate breaking 
of waves rolling in on a gently sloping bottom. Assuming that the 
waves are solitary, we obtain from eq 148 the following limiting 
condition involving the volume Q: 

-^-=>=0 73 

or 

ff=0.7lV5: (215) 

This gives the depth at which the waves will break as a function of 
the volume. If the waves travel up a channel vyith a sloping bed, the 
larger waves break first, and the depths at which they break is pro- 
portional to the square root of their volume. 



The authors express their appreciation to Herbert N. Eaton of the 
National Bureau of Standards for his valuable criticisms and review 
of the paper. 
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